Overview of the Structural, Electronic and Optical Properties of the Cubic and Tetragonal Phases of PbTiO3 by Applying Hubbard Potential Correction

We have performed a systematic study resulting in detailed information on the structural, electronic and optical properties of the cubic (Pm3¯m) and tetragonal (P4mm) phases of PbTiO3 applying the GGA/PBE approximation with and without the Hubbard U potential correction. Through the variation in Hubbard potential values, we establish band gap predictions for the tetragonal phase of PbTiO3 that are in rather good agreement with experimental data. Furthermore, the bond lengths for both phases of PbTiO3 were assessed with experimental measurements, confirming the validity of our model, while chemical bond analysis highlights the covalent nature of the Ti–O and Pb–O bonds. In addition, the study of the optical properties of the two phases of PbTiO3, by applying Hubbard’ U potential, corrects the systematic inaccuracy of the GGA approximation, as well as validating the electronic analysis and offering excellent concordance with the experimental results. Therefore, our results underline that the GGA/PBE approximation with the Hubbard U potential correction could be an effective method for obtaining reliable band gap predictions with moderate computational cost. Therefore, these findings will enable theorists to make use of the precise values of these two phases’ gap energies to enhance PbTiO3’s performance for new applications.


Introduction
Perovskite-type ABO 3 (A = Ba, Sr, Pb) oxides' crystal structure as well as their physical and optical characteristics have been thoroughly explored recently, owing to their widespread development and their functional properties suitable for many technological applications [1][2][3]. In particular, PbTiO 3 itself has been frequently and intensively studied due to its excellent physical properties, which makes it a promising material for a variety of applications, including ferroelectric memory cells [4,5], optoelectronic devices [6,7], dielectric capacitors [8], electromechanical devices [9,10] and pyroelectric applications [11].
Up to now, PbTiO 3 's cubic (Pm3m) and tetragonal (P4mm) phases have been among the most experimentally investigated structures [12][13][14][15]. Indeed, perovskite PbTiO 3 exhibits a ferroelectric-paraelectric phase transition at 763 K. For this reason, a study of the temperature-dependent crystal structure is needed to understand the precise mechanism of ferroelectricity in this material [16,17]. Ferroelectricity is created when cations and anions move away from their equilibrium locations, creating spontaneous dipole moments. The tetragonal phase possesses a non-centrosymmetric space group P4mm and is polar and ferroelectric at low temperatures [18,19]. The spontaneous polarization diminishes as

Computational Model
In the present calculations, we employed the Cambridge Serial Total Energy Package module (CASTEP) [32], which is based on density functional theory (DFT) with the generalized gradient approximation of Perdew-Burke-Ernzerhof [33,34] (GGA-PBE) for exchange-correlation potential. To analyze the interactions involving the valence electrons and ionic core, the valence electron configurations of [Xe] 6s 2 4f 14 5d 10 6p 2 for Pb, [Ar] 3d 2 4s 2 for Ti and [He] 2s 2 2p 4 for O were used. Additionally, self-consistent energy convergence thresholds of 10 −5 eV/atom for total energy, 0.03 eV/ for maximum force, 0.05 GPa for maximum stress and 0.001 for maximum displacement were used to optimize the geometry.
Furthermore, we have systematically examined the dielectric function of PbTiO 3 , which is defined as follows: ε(ω) = ε 1 (ω) + iε 2 (ω) (1) where ω is the incident photon frequency, and ε 1 (ω) and ε 2 (ω) denote the real and imaginary parts of the dielectric function, respectively. To calculate ε 2 (ω), all transitions from occupied to unoccupied states are summed by the following relation [35]: In which e, ε 0 , Ω and ω designate the electronic charge, the vacuum dielectric constant, the volume and the energy of the incident phonon, respectively. Ψ c k and Ψ v k are the wave functions of the conduction band (CB) and the valence band (VB) at the k-point, respectively. The vector that indicates how the incident electric field is polarized is calledû. The momentum operator isû.r and the subscripts c and v denote the CB and VB, respectively [35]. The Kramers-Kronig relation [35] could be used to determine the values of ε 1 (ω): where P is the integral's primary value. Hence, the optical properties using the GGA/PBE + U approximation without the scissor correction, including the optical absorption coefficient (α(ω)), reflectivity (R(ω)), refractive index (n(ω)), extinction coefficient (k(ω)), as well as the energy loss function (L(ω)) were calculated through the complex dielectric function (ε(ω)) and are formulated as follows [35]: where C 0 is the vacuum-state speed of light.
In the geometry optimization section, we will first identify the most appropriate plane wave pseudopotential methods, k-points and cut-off energies to obtain the equilibrium lattice parameters for both the cubic and tetragonal phases of PbTiO 3 with minimal deviations from the experiment via the GGA/PBE approximation. Then, we will search through the Hubbard U potential for the parameters that reproduce the gap energies in perfect agreement with the experimental results. methods for the cubic (Pm3m) and tetragonal (P4mm) phases of PbTiO 3 . While the detailed calculation of the lattice parameters, c/a tetragonality and volume deviations are presented in the Supplementary Materials Section (Table S1).

Figure 1.
Calculation and comparison with the experimental results as a function of the variation in the pseudopotential methods of (a) the lattice parameters and volume of the cubic phase of PbTiO 3 , (b) the lattice parameters (a, c) of the tetragonal phase of PbTiO 3 ; the k-point and cut-off energy values were fixed at 2 × 2 × 2 and 500 eV, respectively, using the GGA/PBE approximation. For both the cubic and tetragonal phases of PbTiO 3 , the pale colors of the lattice parameters and volume indicate that the norm-preserving pseudopotential method has the closest relative deviation to the experimental value. [a] Experimental data from Yoshiasa et al [12]; [b] Experimental data from Shirane et al. [13].
To begin with, we will investigate the appropriate plane wave pseudopotential methods (the ultrasoft, the OTGF ultrasoft, the OTGF norm-conserving and the norm-conserving) for the cubic (Pm3m) and tetragonal (P4mm) phases of PbTiO 3 . For the cubic phase, when we compare with the experimental data [12], the relative deviations of lattice parameters between the calculated structure parameters and the experimental values are about 0.060% for the ultrasoft, 0.075% for the OTGF ultrasoft, 0.256% for the OTGF norm-conserving and 0.055% for the norm-conserving pseudopotential methods. In the case of the tetragonal phase (P4mm) of PbTiO 3 , the results presented in Figure 1b show that the optimized lattice parameters (a, c) provide minimal deviation from the experiment [13] and this is achieved by the norm-conserving pseudopotential method. Hence, the lattice parameters for both the cubic and tetragonal phases of PbTiO 3 calculated by the norm-conserving pseudopotential method display the closest relative deviation to the experimental value.

Appropriate k-Points and Cut-Off Energy
In an attempt to obtain the best minimum deviation from the experimental and optimized lattice parameters, c/a tetragonality, and volume of the cubic (Pm3m) and tetragonal (P4mm) phases of PbTiO 3 , we will investigate the most convenient k-point and cut-off energy values. The pristine cell parameters a = b = c = 3.9702 Å (α = 90 • ; volume = 62.58 Å) [12] and a = b = 3.904 Å, c = 4.152 Å (α = 90 • ; volume = 63.28 Å; tetragonality c/a = 1.063) [13] were originally used for constructing the cubic (Pm3m) and tetragonal (P4mm) phases of PbTiO 3 , respectively. Initially, by keeping the energy cut-off value constant (e.g., 500 eV), we will change the grid values (k-points). Then, after getting an appropriate value of k-points, we will maintain this value constant and vary the cut-off energy to get the most appropriate cut-off energy value. Additionally, further calculation of lattice parameters and volume deviations as a function of varying k-point and cut-off energy values is presented in Tables S2-S5 in the Supplementary Materials Section. In general, the k-points in electronic structure theory are sampling points in the first Brillouin zone of the material, and this set of discrete points in the unit cell allows defining energy values such as energy bands. Setting the k points to 2 × 2 × 2 results in the smallest difference between the experimental and optimized volume and lattice parameters for the cubic phase (Figure 2a coupled with Table S2). Furthermore, following the use of several cut-off energies, we have found that the minimum deviation of the lattice parameters and volume is −0.005% and −0.011%, respectively, occurring for the cut-off energy of 560 eV (Figure 2b coupled with Table S4). On the other hand, for the tetragonal phase, the smallest divergence between the experimental and optimized lattice parameters a and c, as well as c/a tetragonality, occurs by setting the k-points to 2 × 2 × 2 (Figure 2c coupled with  Table S3). Furthermore, after using several cut-off energies, we found that the minimum deviation of a and c for the cut-off energy of 480 eV (Figure 2d coupled with Table S5) is less than 0.65%. Thus, there is good agreement between the optimized tetragonality c/a and the experimental results. For the cubic phase, the distinct colors of the lattice parameters and volume indicate that the k-point and cut-off energy values with the closest relative deviation from the experimental data are 2 × 2 × 2 and 560 eV, respectively. In the case of the tetragonal phase, the pale colors of the lattice parameters (a,c) indicate that the values of the k-points and the cut-off energy that have the closest relative deviation from the experimental results are 2 × 2 × 2 and 480 eV, respectively. [a] Experimental data from Yoshiasa et al [12]; [b] Experimental data from Shirane et al. [13]. Table 1 summarizes the most suitable k-point and cut-off energy values of the cubic (Pm3m) and tetragonal (P4mm) phases of PbTiO 3 using the GGA/PBE approximation. These outcomes highlighted that the difference in lattice parameters between computed and conventional data [12,13] was less than 0.0051% and 0.65% for the cubic (Pm3m) and tetragonal (P4mm) phases, respectively, which allows us to investigate the performance and diver-sity of structural parameters of PbTiO 3 . Hence, the cut-off and k-point values selected after the optimization are in good agreement and demonstrate the validity of our model, together with the smallest relative deviation from the experimental data. Furthermore, in Table 2, a more extensive comparison between our optimized lattice parameters (a, c and c/a tetragonality) and those obtained experimentally [12,13,36,37] and theoretically [25,26,28,[38][39][40], for both the cubic and tetragonal phases of PbTiO 3 is presented. After the geometric optimization of the cubic (tetragonal) phase of PbTiO 3 , we obtained a minimal deviation between the optimized and experimental lattice parameters, when using the norm-conserving pseudopotential method with a 2 × 2 × 2 (2 × 2 × 2) k-point value and a cut-off energy of 560 eV (480 eV). The next step is to perform a good prediction of the electronic structure of the two phases of PbTiO 3 using the Hubbard U correction. Table 1. Computation of the cell parameters, c/a tetragonality, volume and deviations (D) using the GGA/PBE approximation of both the cubic and tetragonal phases of PbTiO 3 with appropriate k-point and cut-off energy values.

Hubbard Potential Correction
For both the cubic and tetragonal phases of PbTiO 3 , it should be noted that the Hubbard U potential of the Pb-5d orbital does not affect the variation in gap energy, owing to its position in the internal electronic structure close to the VB. Furthermore, we demonstrate in the Supplementary Materials Section (Tables S6 and S7) that the gap energy is not influenced even when varying the Hubbard potential values for Pb-5d electrons. To precisely characterize the electronic structures, we used the modified Hubbard U-potential to the Ti-3d and O-2p electrons. The band gaps for the cubic and tetragonal phases of PbTiO 3 using simply the GGA/PBE approximation in the spin polarization domain and without the Hubbard correction are 2.207 eV and 2.213 eV, respectively (Table 3). This difference in band gap shows that the ions (Ti and O) under consideration have strongly correlated electrons, which in turn necessitates that we investigate the Hubbard U correction in our calculations. Table 3. Influence of the Hubbard U parameter values chosen for the O-2p, Ti-3d and Pb-5d orbitals on the gap energies of PbTiO 3 using the GGA/PBE approximation while observing the minimum deviations (D). The cut-off energy/k-point values are 560 eV/2 × 2 × 2 and 480 eV/2 × 2 × 2 for the cubic and tetragonal phases, respectively.
Exp Data. Following the probability series related to the choice of the Hubbard U potential for the Ti-3d and O-2p electrons (see Tables S6 and S7), we have presented in Table 3 and Figure 3 the influence of the chosen Hubbard U parameter (i.e., with minimal deviations) on the band gaps obtained from the tetragonal and cubic phases of PbTiO 3 . Experimentally, the tetragonal phase of PbTiO 3 displays a band gap of 3.40 eV at the Γ-point [14,15], while the band gap of the cubic phase of PbTiO 3 is unknown experimentally. As we can see from Table 3, we fixed the Hubbard U parameters for the O-2p and Ti-3d electrons at 5 eV and 7 eV, and 3.5 eV and 6 eV, for the tetragonal and cubic phases, respectively. This resulted in a gap energy that is exactly 3.400 eV, in conformity with the experimental findings. Moreover, our results also allow us to estimate most of the calculated gap energies which are close to the energies obtained experimentally for the tetragonal phase of PbTiO 3 .

GGA/PBE
In what follows, we will use the Hubbard parameters U p (O) = 5 eV, U d (Ti) = 7 eV and U d (Pb) = 0 eV for the tetragonal phase and U p (O) = 3.5 eV, U d (Ti) = 6 eV and U d (Pb) = 0 eV for the cubic phase to analyze the chemical bonds of PbTiO 3 , including bond lengths, Mulliken charges, effective valence charges and electron charge density. Then, we will also perform a comparative study with and without the Hubbard U correction to investigate the electronic and optical properties of the two phases of PbTiO 3 . experimental findings. Moreover, our results also allow us to estimate most of the calculated gap energies which are close to the energies obtained experimentally for the tetragonal phase of PbTiO3. In what follows, we will use the Hubbard parameters Up(O) = 5 eV, Ud(Ti) = 7 eV and Ud(Pb) = 0 eV for the tetragonal phase and Up(O) = 3.5 eV, Ud(Ti) = 6 eV and Ud(Pb) = 0 eV for the cubic phase to analyze the chemical bonds of PbTiO3, including bond lengths, Mulliken charges, effective valence charges and electron charge density. Then, we will also perform a comparative study with and without the Hubbard U correction to investigate the electronic and optical properties of the two phases of PbTiO3.

Bond Lengths, Mulliken Charges and Effective Valence Charges
The cubic and tetragonal phases of the perovskite-like substance PbTiO3 (with five atoms in the primitive cell) were the focus of this study. In the cubic (Pm3m) phase with local symmetry , the titanium is positioned at the origin of a cubic unit cell, the lead is positioned at the corners of a cube outside the oxygen octahedron, and the oxygen is centered on the face of a cubic unit cell (Figure 4a). Meanwhile, in the tetragonal phase (P4mm) local symmetry , the titanium atom is located in the center of the unit cell, surrounded by six oxygen atoms, and the lead atoms share each of the corners (Figure 4d).

Bond Lengths, Mulliken Charges and Effective Valence Charges
The cubic and tetragonal phases of the perovskite-like substance PbTiO 3 (with five atoms in the primitive cell) were the focus of this study. In the cubic (Pm3m) phase with local symmetry O 1 h , the titanium is positioned at the origin of a cubic unit cell, the lead is positioned at the corners of a cube outside the oxygen octahedron, and the oxygen is centered on the face of a cubic unit cell (Figure 4a). Meanwhile, in the tetragonal phase (P4mm) local symmetry C 1 4v , the titanium atom is located in the center of the unit cell, surrounded by six oxygen atoms, and the lead atoms share each of the corners (Figure 4d).
To determine the bonding nature of Pb-O and Ti-O and to validate our theoretical results based on experimental data, the calculated bond lengths, Mulliken charges and effective valence charges of the cubic and tetragonal phases of PbTiO 3 with GGA/PBE + U approximation are presented in Figure 4b,c,e,f coupled with  [13], as well as with the experimental data obtained by Grazer et al., which are 2.81 Å and 1.99 Å for the Pb-O and Ti-O bond distances, respectively [36], our results are in good agreement with these experimental results as well as with the theoretical calculation (e.g., Ti-O = 1.990 Å) [26]. Furthermore, for the tetragonal phase, the bond lengths of Pb-O1, Pb-O2, Ti-O1 and Ti-O2, are 2.876 Å, 2.557 Å, 1.955 Å and 2.021 Å, respectively (Figure 4e), which agrees well with the experimental data [13].    The Mulliken charge in such a crystal lattice describes the amount of sheared electron density of an atom, with greater positive values indicating that the atom in question contributes greater numbers of electrons; the effective valence represents the difference with the formal ionic charge. The growing covalency shows that the values remain away from zero and perfectly ionic bonding indicates that the effective valence value equals 0. As we can see from Table 4

Electronic Properties
Band Structure and Density of States Figure 6 illustrates the computed band structures in high-symmetry directions in the

Electronic Properties
Band Structure and Density of States Figure 6 illustrates the computed band structures in high-symmetry directions in the Brillouin Zone (BZ) for the cubic and tetragonal phases of PbTiO 3 using the GGA/PBE approximation with and without Hubbard U correction. By applying the GGA/PBE approximation only (Figure 6a,c), an indirect band gap of 2.207 eV and 2.213 eV at the M-G and X-G points is shown as characteristic of the band structure of the cubic and tetragonal phases of PbTiO3, respectively. However, this obtained gap energy value is still much lower than the experimentally obtained one (3.40 eV). Thus, for overcoming the limitations of the GGA, Hubbard U correction was added to the GGA/PBE approximation, as discussed before, resulting in a band gap of 3.400 eV (Figure 6b,d), which is in quite good agreement with the experimental results [14,15]. Because of this, the accuracy of our calculations following optimization and the selection of Hubbard parameters for the cubic and tetragonal phases supports the applicability of the used model.
We computed the total (TDOS) and partial (PDOS) densities of states along with a diagram of the electronic orbital distributions of the cubic and tetragonal phases of PbTiO3 in order to understand the electronic properties, and in particular the atomic orbital contributions in the development of each energy band (Figure 7). With and without the Hubbard correction, the PDOS analysis shows the same electronic contributions (Figure  7a-d) where the VB consists predominantly of contributions from the 2p and 2s orbitals of oxygen and the 5d orbitals of lead, with minor contributions from the 3d and 6s orbitals of titanium and lead, respectively. The 3d titanium and 6p lead orbitals make up the majority of the CB, whereas the 2p oxygen orbitals make up a small portion of it. In addition, the 3p and 4s titanium orbitals contribute to the conduction bands at higher energies. Figure 8 illustrates this graphically, showing an indirect band gap at the M-G Figure 6. Band structure of PbTiO 3 without and with Hubbard potential correction, (a,b) for the cubic phase, (c,d) for the tetragonal phase. Each band energy is shifted with reference to the Fermi level, being set to zero.
By applying the GGA/PBE approximation only (Figure 6a,c), an indirect band gap of 2.207 eV and 2.213 eV at the M-G and X-G points is shown as characteristic of the band structure of the cubic and tetragonal phases of PbTiO 3 , respectively. However, this obtained gap energy value is still much lower than the experimentally obtained one (3.40 eV). Thus, for overcoming the limitations of the GGA, Hubbard U correction was added to the GGA/PBE approximation, as discussed before, resulting in a band gap of 3.400 eV (Figure 6b,d), which is in quite good agreement with the experimental results [14,15]. Because of this, the accuracy of our calculations following optimization and the selection of Hubbard parameters for the cubic and tetragonal phases supports the applicability of the used model.
We computed the total (TDOS) and partial (PDOS) densities of states along with a diagram of the electronic orbital distributions of the cubic and tetragonal phases of PbTiO 3 in order to understand the electronic properties, and in particular the atomic orbital contributions in the development of each energy band (Figure 7). With and without the Hubbard correction, the PDOS analysis shows the same electronic contributions (Figure 7a-d) where the VB consists predominantly of contributions from the 2p and 2s orbitals of oxygen and the 5d orbitals of lead, with minor contributions from the 3d and 6s orbitals of titanium and lead, respectively. The 3d titanium and 6p lead orbitals make up the majority of the CB, whereas the 2p oxygen orbitals make up a small portion of it. In addition, the 3p and 4s titanium orbitals contribute to the conduction bands at higher energies. Figure 8 illustrates this graphically, showing an indirect band gap at the M-G and X-Z points for the cubic and tetragonal phases, respectively, in which the top of the VBs is created by the O 2p orbitals as the bottom of the CBs is formed by the Ti 3d orbitals. From Figure 7, it is evident that the application of the Hubbard U potential correction leads to an increase in the intensity of the 3d orbital of titanium, accompanied by a displacement of 1.187 eV of the CB toward higher energies, which allows restoring the experimental band gap of the tetragonal phase (3.40 eV) [14,15] and this compensates for the systematic error of the GGA approximation. From Figure 7, it is evident that the application of the Hubbard U potential correction leads to an increase in the intensity of the 3d orbital of titanium, accompanied by a displacement of 1.187 eV of the CB toward higher energies, which allows restoring the experimental band gap of the tetragonal phase (3.40 eV) [14,15] and this compensates for the systematic error of the GGA approximation.

Optical Properties
The optical properties of materials have demonstrated themselves to be a valuable means of understanding and confirming the results obtained from the investigation of electronic properties. The following section systematically compares the optical properties of the cubic and tetragonal phases of PbTiO 3 , with and without the Hubbard U correction, with the existing experimental findings. Figure 9 shows the imaginary and real parts of the dielectric function of PbTiO 3 with and without the Hubbard U correction. We will start by discussing the imaginary portion of the dielectric function, which is tightly related to the electronic band structure and shows the light absorption in the material. Using only the GGA/PBE approximation, the ε 2 curve for the cubic and tetragonal phases indicates an energy of~2.207 eV and~2.213 eV, respectively, at the B 0 point (see inset of Figure 9a,c), which is representative of the gap energy of PbTiO 3 . The Hubbard U potential correction to the cubic and tetragonal phases results in a shift of the CB toward higher energies of 1.193 eV and 1.187 eV, respectively. This reproduces the experimental band gap of 3.40 eV at the B 1 point [14,15] (see inset of Figure 9b,d) and corrects the systematic discrepancy of the GGA approximation. From Figure 9a-d, the critical peaks of the optical ε 2 spectra are denoted by A 1 , A 2 , A 3 , A 4 and A 5 , all of which are shifted after the application of the U-correction to higher energies. The transition from O-2p (VBM) to Ti-3d (CBM) orbitals is mostly represented by the A 1 and A 2 peaks, whereas the transition from O-2p (VBM) to Pb-6p (CB) orbitals is primarily represented by the A 3 peak. Moreover, the optical peaks A 4 and A 5 are associated with internal electronic excitation transitions of the Pb-5d and O-2s orbitals near VB to the CB semi-core states. It is necessary to note that each peak in ε 2 is not always connected to a single interband transition because the band structure might contain many direct or indirect transitions with the same energy corresponding to the same peak. Our obtained ε 2 data computed by GGA/PBE + U in the broad energy ranges agree with the currently available experimental findings of the cubic and tetragonal phases of PbTiO 3 in the 0-6 eV range [41,42]; the other transitions were determined on the basis of our obtained results of densities of states ( Figure 7) and agree similarly with the experimental and theoretical data of BaTiO 3 and SrTiO 3 , which have been determined in broad energy ranges [43][44][45][46]. With respect to the real part of the dielectric function, similar results were observed to the imaginary part, where the calculated ε 2 (0) for the cubic (tetragonal) phase are about 6.30 eV (6.53 eV) and 6.95 eV (7.30 eV), respectively, in the case of GGA/PBE and GGA/PBE + U.

Optical Properties
The optical properties of materials have demonstrated themselves to be a valuable means of understanding and confirming the results obtained from the investigation of ranges agree with the currently available experimental findings of the cubic and tetragonal phases of PbTiO3 in the 0-6 eV range [41,42]; the other transitions were determined on the basis of our obtained results of densities of states ( Figure 7) and agree similarly with the experimental and theoretical data of BaTiO3 and SrTiO3, which have been determined in broad energy ranges [43][44][45][46]. With respect to the real part of the dielectric function, similar results were observed to the imaginary part, where the calculated (0) for the cubic (tetragonal) phase are about 6.30 eV (6.53 eV) and 6.95 eV (7.30 eV), respectively, in the case of GGA/PBE and GGA/PBE+U. With and without the Hubbard U correction, the calculated results for the energy dependence of the extinction coefficient k( ) , energy loss function ( ), optical absorption coefficient α(ω) and reflectivity R( ) of the cubic and tetragonal phases of PbTiO3 are shown in Figure 10a,c. As we can see, after the addition of the Hubbard U correction to the GGA approximation, a shift in the optical property spectrum toward higher energy is observed, which leads to correcting and increasing the gap energy of With and without the Hubbard U correction, the calculated results for the energy dependence of the extinction coefficient k(ω), energy loss function L(ω), optical absorption coefficient α(ω) and reflectivity R(ω) of the cubic and tetragonal phases of PbTiO 3 are shown in Figure 10a,c. As we can see, after the addition of the Hubbard U correction to the GGA approximation, a shift in the optical property spectrum toward higher energy is observed, which leads to correcting and increasing the gap energy of PbTiO 3 . From Figure 10b,d, the refractive index n(ω) shows the same behavior after applying the Hubbard U correction. PbTiO3. From Figure 10b,d, the refractive index n( ) shows the same behavior after applying the Hubbard U correction. With respect to the wavelengths (see inset of Figure 10b,d), in order to establish a comparison between our computations and the experimental data, the refractive index of the cubic (tetragonal) phase is 3.27 (2.91) using only the GGA/PBE approximation, and becomes 2.65 (2.41) after applying the Hubbard U correction at the 633 nm wavelength. These results obtained with the Hubbard U correction agree with available experimental findings of refractive indices for sol-gel PbTiO3 films [15] and for metalorganic chemical vapor deposition (MOCVD) PbTiO3 films [6,47,48] by different techniques (Table 5). Based on an earlier examination of the optical characteristics, the Hubbard U potential adjustment was used, which, in turn, effectively offset the systematic error of the GGA approximation and produced results that are in accordance with the experimental data. On the other hand, the electronic analysis carried out for the cubic and tetragonal phases of PbTiO3 is supported by our optical characteristics. With respect to the wavelengths (see inset of Figure 10b,d), in order to establish a comparison between our computations and the experimental data, the refractive index of the cubic (tetragonal) phase is 3.27 (2.91) using only the GGA/PBE approximation, and becomes 2.65 (2.41) after applying the Hubbard U correction at the 633 nm wavelength. These results obtained with the Hubbard U correction agree with available experimental findings of refractive indices for sol-gel PbTiO 3 films [15] and for metalorganic chemical vapor deposition (MOCVD) PbTiO 3 films [6,47,48] by different techniques (Table 5). Based on an earlier examination of the optical characteristics, the Hubbard U potential adjustment was used, which, in turn, effectively offset the systematic error of the GGA approximation and produced results that are in accordance with the experimental data. On the other hand, the electronic analysis carried out for the cubic and tetragonal phases of PbTiO 3 is supported by our optical characteristics.

Conclusions
In conclusion, this work presents a comprehensive overview of the crystal structure, chemical bonds, electronic and optical properties of the cubic and tetragonal phases of PbTiO 3 using the GGA/PBE approximation with and without Hubbard U correction. First of all, we performed detailed calculations regarding geometry optimization to obtain the smallest difference across the experimental and optimized lattice parameters/volume, with the selection of the most appropriate methods of plane wave pseudopotential, kpoints and cut-off energy, which validates the accuracy of our model. Secondly, the electronic investigations have shown that adding the Hubbard correction to the GGA/PBE approximation, with the well-selected values of Hubbard parameters for p-and d-orbitals, achieves an excellent result for the band gap (3.400 eV) of the cubic and tetragonal phases of PbTiO 3 as compared to the experimental data. Additionally, the optimized chemical bonding analysis shows consistency with the experimental results, suggesting that the Pb-O and Ti-O bonds have a partial and strong covalent bonding nature, respectively. In addition, the application of the Hubbard U potential correction to the calculated optical properties could overcome the systematic inaccuracy of the GGA approximation. At the same time, studies of the complex dielectric function and other optical properties confirm the electronic calculations and offer excellent concordance with the actual results over a wide range of energies. Consequently, these results can be leveraged by scientists to expand the targeted scope of the cubic and tetragonal phases of PbTiO 3 .

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ma16124302/s1, Table S1: Computation of the cell parameters, c/a tetragonality, volume and deviations using the GGA/PBE approximation of the cubic (Pm3m) and tetragonal (P4mm) phases of PbTiO 3 as a function of the variation in the pseudopotential methods; the k-point and cut-off energy values were fixed at 2 × 2 × 2 and 500 eV, respectively; Table S2: Computation of the cell parameters and volume deviation using the GGA/PBE approximation of the cubic phase (Pm3m) of PbTiO 3 as a function of the variation in the k-point values; the value of cut-off energy was set at 500 eV; Table S3: Computation of the cell parameters (a, c) and tetragonality c/a deviation using the GGA/PBE approximation of the tetragonal phase (P4mm) of PbTiO 3 as a function of the variation in the k-point values; the value of cut-off energy was set at 500 eV; Table S4: Computation of the cell parameters and volume deviation using the GGA/PBE approximation of the cubic phase (Pm3m) of PbTiO 3 as a function of the variation in the cut-off energy values; the k-point values were set at 2 × 2 × 2; Table S5: Computation of the cell parameters (a, c) and tetragonality c/a deviation using the GGA/PBE approximation of the tetragonal phase (P4mm) of PbTiO 3 as a function of the variation in the cut-off energy values; the k-point values were set at 2 × 2 × 2; Table S6: Hubbard U parameter values chosen for the O-2p, Ti-3d and Pb-5d orbitals of the cubic phase (Pm3m) of PbTiO 3 using the GGA/PBE approximation; Table S7: Hubbard U parameter values chosen for the O-2p, Ti-3d and Pb-5d orbitals of the tetragonal phase (P4mm) of PbTiO 3 using the GGA/PBE approximation.
Author Contributions: I.D., performing theoretical calculations, methodology and writing-original draft preparation; M.A., validation and writing-review and editing; R.I.E., validation and writing-review and editing; A.I.P., validation and writing-review and editing; A.R., validation and supervision. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
The data presented in this study are available from the corresponding author upon reasonable request.